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We compute the coefficients of bulk viscosity for a non-relativistic superfluid corresponding to a 
fermionic system close to the unitarity limit. We consider the low temperature regime assuming 
that the transport properties of the system are dominated by phonons. To compute the coefficients 
of bulk viscosity we use kinetic theory in the relaxation time approximation and the low energy 
effective field theory of the corresponding system. We show that the three independent bulk viscosity 
coefficients, (^i,("2,C3i associated with irreversible flows vanish for phonons with a linear dispersion 
law. Considering a phonon dispersion law with a cubic term in momentum we find that in the 
conformal limit Ci = C2 = 0, while ^3 is non-zero. Including a conformal breaking term which arises 
for a large but finite s-wave scattering length, a, at the leading order in 1/a we obtain that i^i oc 1/a 
and ("2 oc 
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I. INTRODUCTION 

The properties of quantum degenerate fermionic systems with an attractive two-body scattering interaction have 
been the subject of extensive investigation in the last years Of particular interest are systems with an infinite 
two-body scattering length that are believed to have universal properties [2| , meaning that the features of the system 
are independent of the detailed form of the inter-particle potential. 

Experiments with trapped cold atomic gases are able to reach the region of infinite scattering length (the so-called 
unitarity limit) tuning the interaction between the fermionic atoms by means of a magnetic-field Feshbach resonance 
[1) 0) H) @]- In general, in these experiments the two different populations of fermions consist of atoms, like ^Li or 
^''K, in two hyperfine states. The strength of the interaction between atoms depends on the applied magnetic field 
and can be measured in terms of the s-wave scattering length. By varying the magnetic-field controlled interaction, 
fermionic pairing is observed to undergo the Bose-Einstein condensate (BEC) to Bardeen-Cooper-Schrieffer (BCS) 
crossover. The unitary limit is reached when the magnetic field is tuned at the Feshbach resonance 0, where the 
two-body scattering length diverges. 

Far from the unitarity limit the properties of the system are qualitatively well understood using mean field theory . 
In the weak coupling BCS region the system is characterized by the formation of Cooper pairs. In the strong coupling 
limit the system can be described as a BEC dilute gas. The extreme BCS and BEC regimes are also in good 
quantitative control in mean field theory. However, the mean field expansion is not reliable close to unitarity because 
the scattering length is much larger than the inter-particle distance and there is no small parameter in the Lagrangian 
to expand in. Therefore fluctuations may change the mean field results substantially. 

Close to the unitarity region a quantitative understanding of the phases comes mainly from Monte-Carlo simula- 
tions . Other approaches consist in considering the expansion in a small parameter that comes from the generaliza- 
tion to an arbitrary number N of spins [l^l, or to d space dimensions [ll|. In the former case, for N ^ 00 the problem 
is exactly solvable by mean-field theory and one can consider corrections and then extrapolates to = 2. In 
the second case one considers an e = 4 — d expansion and then extrapolates the result to d = 3. Both approaches are 
in quantitative agreement with Monte-Carlo simulations for the equation of state of the system. A different method 
consists in adding the quantum fluctuations on the top of the mean- field theory [l3,[l3. In this way one improves 
the agreement with Monte-Carlo simulations, with respect to the bare mean-field calculation. 
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It is a remarkable aspect of these fermionic systems that for any value of the attractive interaction they are superfluid, 
provided the temperature is sufficiently low. Superfluidity is a phenomenon that occurs after the appearance of a 
quantum condensate that breaks a global U{1) symmetry of the system [13, [H, [H] , regardless of whether the system 
is fermionic or bosonic. In both cases Goldstone's theorem predicts the existence of low energy modes with a linear 
dispersion law, which are essential to explain the property of superfluidity. We will refer generically to these modes 
as superfluid phonons, or phonons for simplicity, which dominate the transport effects of the system at very low 
temperature. 

The hydrodynamic equations governing the bulk fluctuations of a superfluid arc essentially different from standard 
fluid equations. At non- vanishing temperature one has to employ the two-fluid description of Landau p^ . which takes 
into account the motion of both the superfluid and of the normal component of the system. In order to describe the 
different dissipative processes one has to introduce more transport coefficients than in a normal fluid. In particular, 
one has three independent bulk viscosities Ci: C2) Csi ^is well as the shear viscosity and the thermal conductivity. 
The shear viscosity of a unitary superfluid at low temperature has been computed in Ref. [l^. In Ref. [l^ it has 
been shown by a general argument that (1 and C2 vanish in the conformal limit [l9l |. However, ^3 and the thermal 
conductivity associated with irreversible heat flow cannot be determined by the same symmetry reasoning. 

In the present paper we evaluate the three bulk viscosity coefficients in the low temperature regime, T <C Tc, 
where Tc is the critical temperature for superfluidity, as we consider the contribution of phonons only, assuming that 
the contribution of other degrees of freedom is thermally suppressed. We show that all the bulk viscosities vanish 
for phonons with a linear dispersion law. This is a result that was obtained several years ago by Khalatnikov and 
Chernikova [2^, although not widely known. These authors also found that for superfluid ^He all the transport 
coefficients, with the exception of the shear viscosity, vanish if one considers only phonons with a linear dispersion 
law. Then, we evaluate the transport coefficients considering a phonon dispersion law that includes cubic corrections. 
We report the general expressions for the three dissipative coefficients, given in terms of the parameters that appear in 
the phonon dispersion law and in terms of the pertinent decay rate. Our results strongly depend on the coefficient of 
the cubic term in the phonon dispersion law, which is poorly known. If experiments and/or Monte-Carlo simulations 
could measure more precisely the value of this coefficient then we would know with more accuracy the numerical 
values of the bulk viscosities. 

For the specific system we are interested in, the properties of the phonons needed in our computation can be 
extracted from the effective field theory constructed in Ref. [2l|. The Lagrangian of this theory is determined by 
demanding non-rclativistic general coordinate invariance and conformal invariance and assuming that phonons are 
the only relevant degrees of freedom. From the effective Lagrangian one obtains that the phonon dispersion law and 
self-couplings depend on some universal and dimensionless constants. Employing the expressions for these quantities 
we show that in the conformal limit Ci = C2 = 0, while (^3 is non-zero. 

As a final step we study a system that is close to unitarity. In this case scale invariance is broken and additional 
terms in the effective Lagrangian must be taken into account [2l| . In the presence of a large but finite s-wave scattering 
length a, we study how the phonon dispersion law and the three-body self-couplings are modified. Then, we evaluate 
the first non- vanishing corrections to the bulk viscosity coefficients and find that Ci cx 1/a and (,2 oc 

It is also interesting to compare our results with those corresponding to a Bose superfluid. A computation of all 
the transport coefficients for a dilute but condensed Bose gas due to phonons was presented in Ref. (2^, see also 
Ref. [l6| for a more extended discussion. Remarkably, the temperature dependence for the bulk viscosities in these 
two different superfluids are the same. 

Let us finally point out that the techniques we employ, and even the explicit computations, are very similar to those 
used in the evaluation of the transport coefficients for relativistic superfluids, in particular for the color flavor locked 
phase of dense quark matter (2^, Hjl . 

This paper is organized as follows. In Scction[TT]wc briefly review the superfluid hydrodynamics for a non-relativistic 
ffuid including effects due to dissipation. The dissipative coefficients can be determined using kinetic theory and in 
Sectionlllllwe give general expressions for the bulk viscosity coefficients employing the two different methods described 
in Ref. [1^. In Section [TVl we present the low energy effective theory valid for a cold Fermi gas, in the exact unitarity 
limit. Sec. IIV Al and close to it. Sec. IIVBI The explicit evaluation of the bulk viscosity coefficients is reported in 
Sec. |Vl We draw our conclusions in Sec. IVII 

Throughout, we use natural units, so that we take the Boltzmann and Planck constants as fc^ = ?i = 1 in all our 
computations, but we restore them in the final section. 



II. SUPERFLUID HYDRODYNAMICS AND KINETIC THEORY FOR THE PHONONS 

In this Section we briefly review the two-fluid theory of superfluidity developed by Landau. We also provide the 
main ingredients of the kinetic theory for superfluid systems described by Khalatnikov [l^ . This Section might be 
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skipped by the reader who is famihar with supcrfluid hydrodynamics. 

A. Hydrodynamics of a non-relativistic superfluid 

In a superfluid there are two independent motions, one normal and the other superfluid. with velocities v„ and v^, 
respectively. These motions are associated with two different matter and current densities, such that the total density 
and the total current density of the system are given by the sum of the superfluid and normal components as 

P = Pn+ Ps , j = /3«V„ + PsVs . (1) 

The superfluid motion is irrotational, curlv^ = 0, and thus its velocity can be written as the gradient of a scalar 
function that is proportional to the phase of the wave function condensate. 

The hydrodynamic equations have the form of mass and momentum conservation laws and neglecting dissipation 
they are respectively given by 

dtp + divj = , (2) 

and 

dtU + ^kn^k = , (3) 

where 

n^fe = PnVmVnk + PsVsiVsk + P^ik , (4) 

is the momentum flux density tensor, and P is the pressure of the system. 

Since in a superfluid there can be two different motions, beside Eq. ^ one has a second hydrodynamic equation 
describing the irrotational motion of the superfluid component 

a^v. + V + y) =0, (5) 

which indicates that a gradient in the chemical potential acts as a force for the superfluid component. 
In order to complete the system of equations, the energy conservation law is also needed 

dtE + div(Q) , (6) 

where E is the energy per unit volume and the energy flux is given by 

Q = (^M + j + ST^f.„ + p„v„(v„ • (v„ - V,)) , (7) 

where S is the entropy. In the absence of dissipation, entropy is conserved and one has that dtS + div(5vn) = 0. 

In the presence of dissipative processes there arc additional contributions to the hydrodynamic equations that arise 
from irreversible processes, thus 

^+9,(%+r,,) = 0, (8) 
^ + V.(Q + Q') = 0, (10) 

where 

Q' = q + /i(j -pv„) + r •v„, (11) 

and Ty , h and q are small dissipative terms. In this case entropy is not conserved and the entropy production rate is 
given by 

i?= - W- (p,(v„- V,)) -T,fc5fc?;„,- ^q- VT. (12) 
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From the requirement that the dissipative processes induce an increase of the entropy it follows that 

2 

Tij = -?y(9ji'„.j + a^Wnj - -(5yV • v„) - (5y (CiV • (ps(vs - v„)) + C2V • v„) , (13) 

h = -C3V • (p,(v, - v„)) - C4V • v„ , (14) 
q = -kVT, (15) 

where 77 is the shear viscosity, k is the thermal conductivity associated with irreversible heat flow, and Ci j C2 1 Ca i C4 ^-re 
the four bulk viscosity coefficients. According to the Onsager symmetry principle, the transport coefficients satisfy 
the relation ^1 = C4j while the requirement of positive entropy production imposes that k, rj, (^2, C3 ^-re positive and 
that Cf < (2(3- 

Due to the presence of various bulk viscosity coefficients, the liydrodynamic equations are quite involved. From 
Eg. p^ it is clear that (2 plays the same role as the standard bulk viscosity coefficient. On the other hand, from 
Eqs. (|13p and (jl4p . one can see that Ci, Ca and ^4 provide a coupling between the hydrodynamic equations of the two 
components. It is the presence of these couplings that makes the hydrodynamic equations complicated. 

However, notice that the force acting on the superfluid component is proportional to the gradient of h and therefore, 
even in the presence of dissipation, the motion of the superfluid component will be almost everywhere irrotational, 
with vorticity concentrated in superfluid vortices, see e.g. [25j . 

Let us assume that the superfluid component has a constant density. The dissipative terms associated with Ci and 
(■3 vanish when Vg — = const. Therefore, these viscosity terms oppose to the presence of arbitrary, space dependent, 
relative motion between the two components. However, the motion of the normal component is not irrotational (due 
to the presence of the shear viscosity) whereas the superfluid component is irrotational. This means that it is never 
possible for these dissipative terms to be zero, unless d and ^3 are zero or unless a superfluid vortex is created. 
Finally, the effect of the dissipative term proportional to ^4 is to produce a force on the superffuid component, due to 
the variation of the velocity of the normal component. 

The friction forces due to bulk viscosities can also be understood as drops in the main driving forces acting on the 
normal and superffuid components. These forces are given by the gradients of P and fi, respectively. Therefore we 
can write that 

P = Pcq - Cldiv(ps(Vn - Vg)) - C2divVn , (16) 

^J■ = A'cq ^ Ca div(ps(v„ - Vg)) - C4divv„ , (17) 

where Poq and /icq are the pressure and chemical potential in the absence of bulk viscosities. We shall use this 
interpretation for the computation of the bulk viscosity coefficients. 

In the following Sections we will find that in the conformal limit C4 — Ci ~ C2 = whilst (^3^0, meaning that 
the only dissipative force acting on the normal component is due to the shear viscosity. This force will be a source of 
vorticity for the normal component and since (^3 is non-zero, it will generate a force acting on the superffuid component 
that will tend to make Vg — Vn constant. 

B. Kinetic theory for the superfluid phonons 

At very low temperatures phonons give the leading thermal contribution to all the thermodynamic properties of 
the superffuid. In the hydrodynamic regime, phonons also give the leading contribution to the transport coefficients 
entering into the two-ffuid equations. Khalatnikov developed the kinetic theory associated to these degrees of freedom 
that we brieffy review [T5| . 

Let us first note that at higher temperatures various quasiparticles or collective modes may contribute to the ther- 
modynamics or to transport phenomena. For superfluid ^He the energy spectrum reveals the presence of excitations 



called rotons, that can be taken into account in the construction of the kinetic theory [15[ . For the cold Fermi gas close 
to the unitary limit, finite temperature Quantum Monte-Carlo simulations [2^ and experimental measurements [27j . 
reveal that only at very low temperatures the thermal spectrum might be well-described by phonons. At higher 
temperature the spectrum is much more complicated, and in the computation of transport coefficients one may need 
to include other different contributions. However, we shall not consider such a high temperature regime. 
In general one can assume that the dispersion law of phonons is given by 

ep = c,p + V + 0(p^), (18) 



where Cs is the speed of the phonon and for systems with a small coefficient of thermal expansion it is equal to the 
speed of first sound. 
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Under the assumption that the leading contribution is due to phonons, one can compute various thermodynamic 
quantities starting from the phonon distribution function n. The entropy density is given by 

V = ^ / dppMn + , (19) 



the number of phonons per unit vohime is given by 



A/'ph=/(03^' (20) 



and the phonon contribution to the total chemical potential is given by 

^"^"^i/"^^'"^^' ^^^^ 

where /Uq is the chemical potential at zero temperature. 

In the study of transport phenomena it is necessary to consider the evolution of the out-of-equilibrium phonon 
distribution function n, which obeys the kinetic equation 

dn ^ dn dH dn dH 

dt dr dp dp dr ' 

where H = ep + p ■ Vs is the phonon Hamiltonian, and C[n\ is the collision integral. 
At equilibrium phonons follow the Bose-Einstein distribution 

"cq(ep) = ge,/T _ 1 ' (23) 

and the collision term vanishes. For small departures from equilibrium one can linearize the collision term on the 
deviations 6n = n — Ucq and the transport coefficients can be obtained by solving the kinetic equation obeyed by Sn. 
In general, this is a complicated task, as one has to deal with an integro-differential equation. For our purposes it is 
sufficient to obtain an approximated expression for the transport coefficients and therefore we shall use the relaxation 
time approximation (RTA). In the RTA the collision term is written as 

SC = -^, (24) 

Trcl 

where Trci is the relaxation time for the coUisional process that gives the leading contribution to the transport 
phenomena one is studying; for the bulk viscosity coefficients, collisions that change the phonon number. 

Using the RTA one can easily obtain the solution for Sn, and the corresponding dissipative fluxes in the hydrody- 
namic equations. The RTA provides the correct parametric dependence of the various transport coefficients on the 
relevant scales, although it does not fix with accuracy the numerical factor in front of these quantities. In the present 
article we shall content ourselves with the RTA. The reason is that for a system of cold fermionic atoms the coefficient 
of cubic order in the phonon dispersion law has not been precisely determined, therefore there is little motivation for 
achieving a good precision on the numerical factors of the transport coefficients. 



III. PHONON CONTRIBUTION TO THE BULK VISCOSITY COEFFICIENTS 



In this Section we compute the bulk viscosity coefficients for a non-relativistic superfluid. Khalatnikov described 
two different methods for the evaluation of these quantities, one based on studying the evolution of the phonon number 
density Afph, see Sec. IIII Al and the other one based on studying the evolution of the phonon distribution function 
n, see Sec. IIII Bl Here we present in detail both methods and show that the first method corresponds to solving the 
transport equation in the relaxation time approximation, and thus, it is equivalent to the second method. 

We find that for phonons with a linear dispersion law all the bulk viscosity coefficients vanish, independent of 
whether the system is conformal invariant or not. Then, we evaluate the bulk viscosity coefficients for phonons with 
a cubic dispersion law and express the result in terms of the parameters B and Cs, see Eq. (jl8p . 
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A. Evaluation of the bulk viscosity coefficient with the first method 



When a perturbation applied to a superfluid system determines a change of the number of phonons per unit vohime, 
A/'ph, colhsional processes tend to restore the cquihbrium value of this quantity. The evolution equation for the phonon 
number can be written as 

dtAfpb + div(7VphV„) = — Y^Hph , (25) 

where the rate of change is expressed as a power expansion in a "fake" phonon chemical potential, /iph, and the 
decay rate of phonon changing processes, Fph [28j . Expressing the phonon number as a function of the density and 
of entropy, and using the linearized continuity equations for these quantities, one can express the phonon chemical 
potential in terms of the different dissipative flows that appear in the hydrodynamic equations. These terms modify 
the equilibrium pressure and chemical potential, and comparing the results with the expression in Eqs. p6|) and (jl7p . 
one identifies the different bulk viscosity coefficients. 

For small departures from equilibrium and for small values of v^ and v„ it turns out that [l5l | 



and therefore = C2C31 rneaning that one of the relation for positive entropy production is saturated. Here, we have 
defined the quantity 

A = ^ , (29) 
dp 

while 

h=Mp^-S^-p^, (30) 

and in order to evaluate the various derivatives that appear in these expressions we change variables. Consider that 
in Eqs. ([29|l and ([30|l it is assumed that the independent variables are S and p. Now, we write S ~ S{T,po) and 
p = p{T, pq), and by the chain-rule we have that 

(31) 
(32) 

up ui up upo up 

One can simplify these expressions 
using the Maxwell relation 

(33) 

and it is now easy to check that with a linear dispersion law all bulk viscosity coefficients vanish. In order to evaluate 
the leading correction in B to the viscosity coefficients we define the adimensional parameter 

X = ^ , (34) 

and expand the various quantities evaluated with the equilibrium phonon distribution function to the leading order 
in X. In this way we obtain that the number of phonons per unit volume is given by 

- 2^ (r(3)C(3) - X r(6)C(5) + Oix')) , (35) 
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while the entropy turns out to be 

^f" = ^ (r(5)C(4) - 3xr(7)C(6) + Oix^)) , (36) 
and the chemical potential is given by 

" ^ + aSj (t^-CXC) + =='-(fiK(6)(|f - 6f ) + o(.')) , (37) 

In all the above expressions T(z) and C(z) stand for the Gamma and Riemann zeta functions, respectively. 
From these expressions we have that the first non-vanishing corrections to Ii and I2 are 

^.^^-=('^«3,-7«.,)(e.f-.|.), 

and 

h = -^T^(-^C(3) - 7C(5)) - pi, = -y^T^(-^C(3) - 7C(5)) (2BC. + 3p (c.^ - B^)) ■ ^ 

These expressions are fully general, valid for any non-relativistic superfluid to the leading order in x. Notice that 
for the next-to-leading order temperature corrections one has to include terms in the phonon dispersion law going as 

, that we have neglected. In this paper we will only consider the first non- vanishing correction to the bulk viscosity 
coefficients, the next-to-leading order corrections being very suppressed at low temperatures. 

Once the explicit expressions of Cs and B and their dependence on the density are known, Ii and I2 can be evaluated. 
The quantity that remains to be evaluated is the decay rate of phonon changing processes, Fph. The parameter B 
determines whether some processes are or are not kinematically allowed. For B > the leading contribution comes 
from the Beliaev process </> — *■ 00. In the opposite case one has to consider processes like (pcj) — > 00</> [2§| . 



B. Evaluation of the bulk viscosity coefficients with the second method 

The transport coefficients that appear in the superfluid equations can be derived by studying the evolution of the 
deviations from equilibrium of the phonon distribution function, Sn. The Boltzmann equation (|22p . amended with a 
collision term describing phonon number changing processes, is linearized to get the equation obeyed by 6n. At this 
stage one could use, for example, the Chapman-Enskog procedure, see e.g. [221 ■ This consists in assuming that dn 
can be expressed as a function of the hydrodynamical variables and their grandients. Then, the equation could be 
solved by a variety of numerical approaches (variational method, use of orthonormal polynomials, etc). At a technical 
level, the problem is very much simplified with the use of the RTA. Then one gets a simple analytical solution for Sn, 
which is good enough to get an approximated value of the transport coefficients. The RTA allows us to get the correct 
parametric dependence of the transport coefficients on the scales of the problem, but it does not fix with accuracy 
the numerical factor in front of these quantities. 

Here we show that using the RTA we obtain the same results found in the previous Section if we take the relaxation 
time as 



1 Tph 

Trcl A/'ph ' 



(40) 

^Vph 

where Fph is the phonon decay rate appearing in Eq. (|25p . 

With the use of the RTA one finds the following expressions for the bulk viscosity coefficients [l^ 

ci = -r.c. / -'^-^^ J, (}p- ^ , (41) 

G = -r.e: / --^^-" J2 (}P- ^ + . (42) 
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where we have defined the quantities 

, SdT lde„ 
J. - PJ^ + TOs'p--,^-P- (46) 

At the leading order in B all the integrals are zero, so we evaluate Ji and J2 at order B^. Then we realize that 
they can be written as 

/ dB ^dcs\ M 

f f dB dcsW M , , 

where we have defined the quantities 

M = ~lc%v'^ ^YIQ'k^BcIp^T'^ + 2Qc\-K^T'^ -W)%t:^Bc,T^ ^0{B'^) , (49) 
L = lcl~\mBc\T{^T'^ ArO{B^). (50) 

After performing the final integral in momenta, we find that the expressions of the bulk viscosities obtained by the 
two methods have the same parametric dependence on physical quantities. While this is what one could naturally 
expect, this result only comes out after several subtle cancellations of different terms in the expansion. 

Although we are not interested in precise numerical factors, we notice that taking 

(51) 

the two methods give the same numerical expressions for the bulk viscosities. In this way, we check that the first 
method presented here is equivalent to solving the transport equation in the RTA. 

IV. LOW ENERGY EFFECTIVE FIELD THEORY FOR THE GOLDSTONE MODE 

The information required for the evaluation of the bulk viscosity coefficients might be extracted from the effective 
field theory associated to the Goldstone mode of the superfluid system. The effective field theory is constructed as an 
expansion over derivatives and over Goldstone fields, and it provides a systematic way to compute different physical 
quantities to a given accuracy. 

It has been known for a while that the leading order term Lagrangian of the Goldstone mode of a superfluid 
system is entirely fixed by the equation of state [30,[ll|- For the unitarity Fermi gas, the next-to-leading Lagrangian 



has been constructed in Ref. 21[, requiring that it is invariant with respect to non-relativistic general coordinate 
and conformal transformations. We review here the low energy effective field theory for the cold Fermi system, and 
extract from it the parameters required for the evaluation of the bulk viscosity coefficients, both in the exact unitarity 
limit fSec. lIV A|) . and close to the unitarity limit (Sec. IIV B|) . We would like to stress that the methodology used here 
might be also followed to study other superfluid systems, governed by different equations of state, and with different 
global symmetries. 

A. Effective Lagrangian in the exact conformal limit 

In the unitarity limit, the thermodynamic properties of the cold Fermi gas can be determined up to some dimen- 
sionless constants At zero temperature, and due to the absence of any internal scale, dimensional analysis fixes 
the form of the pressure as being proportional to that of a free system 

P = comVo^' , (52) 
where cq is a dimensionless and universal constant. This parameter can be expressed as 

25/2 

^0 = 1^^2^' (53) 
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where ^ is the universal constant that fixes the relation between the chemical potential and the Fermi energy fiQ = 
S^Ep/m. Experiments with cold trapped fermionic atoms [s^j find f_~ 0.32 — 0.44, a result that is in agreement with 
Quantum Monte-Carlo calculations at vanishing temperature [s^ [sl] . 
The density at zero temperature is then easily deduced 

4 3/2 

Po = -^com , (54) 

and the speed of sound at T = turns out to be 




= \ Poir- = \ Hr ■ (55) 



The leading order (LO) Lagrangian for the Goldstonc field can be determined by the equation of state; the reason 
being that the effective action of the theory at its minimum for constant classical field configurations has to be equal 
to the pressure. The next-to-leading order (NLO) Lagrangian is constructed by demanding invariance with respect 
to non-relativistic general coordinate invariance and conformal invariance [2l| . The combined terms then read 

C = C^o + C^LO = com'^'X'/^ + cimV2 IZ^ + ^(v^^)^^^, (56) 

vx v"i 

where 

X = mfio~ ao</. - , (57) 

Zm 

and (j) is the phase of the condensate and we have neglected effects due to the trapping potential. Notice that /io is 
a chemical potential with dimensions of velocity squared. Our definition of chemical potential is the one given by 



Khalatnikov in Ref. [15| and differs from the definition given in Ref. [21| by a mass factor. Notice that at this order 
the effective Lagrangian depends on two more dimensionless parameters ci and C2, that are universal constants. 

If one expands the effective Lagrangian in Eq. (|56p in the field (f) up to quadratic order, and brings the kinetic term 
into standard form, one obtains that the phonon dispersion relation up to cubic powers of momentum is given by (2]| 

ep = Cs (^P-TT^V^ (^ci + ^C2^ , (58) 

where kp is the Fermi momentum, Ep = kp/2m. Then, the speed of the phonon agrees with the speed of sound, 
Eq. ([551) . ^ii'i can identify the value of the parameter B introduced in Eq. as 

B=~TT^C,^(^Cl + ^C2^^. (59) 

The two parameters Ci , C2 are related to the momentum dependence of the static density and transverse response 
functions and can be evaluated by the e-expansi on [35l . finding ci ~ —0.0209, and C2/C1 = 0{e^). These parameters 
can also be evaluated within mean-field theory [30], finding a value for ci that differs from the previous one by a 
30%. While the numerical discrepancy between the predictions of the two methods is large, both give Ci -I- |c2 < 0, 
meaning that B is positive and the Bcliaev process (f> — > (fxf) is kinematically allowed. We will assume that the sign of 
B is correctly predicted by the two methods, but leave ci and C2 as coefficients still to be determined. 

The three-phonon self-coupling is also determined by expanding the Lagrangian above, and reads [l8l | 

£30 = ~a {{docl>f - 9c2(ao0)(V0)2) + . . . , (60) 
where terms with higher number of derivatives have been neglected, and the coupling is given by 



TTC 



3/2^3/4 



(61) 



The effective field theory here presented allows us to compute the corrections to the phonon dispersion law needed in 
our computations. We have checked that the one-loop corrections to the dispersion law are very suppressed. Further, 
they do not modify the sign of the coefficient B that appears at tree level, which is relevant in deciding whether the 
Beliaev decay is kinematically allowed or not. A more detailed discussion of the one- loop phonon dispersion law will 
be presented elsewhere. 
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B. Effective Lagrangian close to the conformal limit 

For finite scattering length, the pressure can be written as a power expansion in 1 j a 

P = Po + PcB = comV^/^ + + • • • , (62) 

a 

where do is a dimensionless constant to be determined by matching or experimentally. 

This change in the equation of state induces new terms in the efi^ective Lagrangian which are responsible for the 
breaking of scale and conformal invariance [2l| 

£ = £lo + C^LO + £cB + • • • = com^/^X^/' + cimV^iZ^ + ^(v^^)^ + do"^ + ■ ■ ■ (63) 



The effect of the scale breaking term is to change the phonon dispersion law, and the self-coupling coefficients. In 
order to find the proper transformation of these quantities we expand £ to quadratic order in 0, and bring the kinetic 
term into the canonical form by the field rescaling 

- - (1 - y) , (64) 



where we have defined the quantity 



y = 7S= ■ (65) 



One then finds that the speed of the phonon is corrected at the leading order in 1/a as 



c.,CB = (l + I) , (66) 
while the cubic term in the phonon dispersion law is now given by 

(-1)4=^-1))- 

The three-phonon self-coupling is also modified as 

£cB = -«CB [{docP)' - 9clcB (1 + y) [dom^^f] +■■■ (68) 

where the coupling reads 

acB = a (1 - 3y) . (69) 
The density at zero temperature of the system also deviates from its value in the exact unitarity limit by 

Po^=Po(l + 3y) . (70) 

These corrections to the parameters of the phonon dispersion law and to the self-couplings coefficients will turn out 
to be of great importance in the computations of the bulk viscosity coefficients, because they will allow us to compute 
the only non- vanishing contributions to Ci and C2 close to the conformal limit. 

V. BULK VISCOSITIES CLOSE TO THE UNITARITY LIMIT 



For the explicit evaluation of the bulk viscosity coefficients, we start by considering the exact unitarity limit. Using 
the expressions of and B given respectively in Eqs. (|55p and (|59p . we evaluate the derivatives with respect to the 
density. At low temperature we can approximate 



11 



where we have used the fact that = |^ . In a similar way wc find that at unitarity and for vanishing temperatures 

Upon substituting the expressions in Eqs. ([7T|) and ([7^ in Eqs. ([55]) and ([5^ it follows that 
and 

/2 = . (74) 

Therefore the bulk viscosity coefficients Ci and ^2 vanish. 

To obtain the final expression for we still have to evaluate the decay rate Fph- If the process (f) (fxf) is 
kinematieally allowed, as suggested by the results of Refs. [s^, [s^, then we have that 



^P'^^ / {2^f {2nY {2^Y I^^P"oq(ep) (1 + neq(e,)) (1 + neq(6,)) {27:)H^^\P K Q) , 



(75) 



where the delta functions ensure the energy- momentum conservation, where P = (ep, p), Q = (e^, q), K = (e^, k) and 
where |M(P, Q, K)\'^ is the squared of the scattering amplitude with non-relativistic normalization of the integration 
measure. This quantity is related to the squared of the scattering matrix with relativistic conventions |A^p by 

|M(P.q,A-)r J^'7f)l\ (70, 



and can be computed from the three-phonon interaction, Eq. ()60|1 . The scattering matrix takes a particular simple 
form if expressed in terms of the vector product defined as P ■ K ~ €pek — 9c^p • k [l^ 

\MiP, Q, K)f = {epQ-K + eqP-K + e,. Q ■ Pf , (77) 

where a is defined in Eq. (j6ip . 

For phonons with a linear dispersion law, the splitting process is perfectly collinear. For B > the process is non- 
coUinear, and the angle of scattering between the incoming phonon and the two outcoming phonons is proportional 
to B. Since Ii and I2 are proportional to B we can evaluate the decay width considering the collinear process only 
and from the expressions above we obtain that 

Fph ~ 14781.6^-^ + 0{x) , (78) 

where x is defined in Eq. (|34p . Note that in the exact unitarity limit x ^ (T//io)^, so that at low temperature it is a 
good approximation to neglect terms of order x in Fph. 

Upon substituting the decay width and the expression for Ii reported in Eq. ([75)1 in Eq. ([^5]) we find that 

^3 0.015^?^ + O (T') ~ 3695.4 f ^ V^' (^i±|^y3 ^ ^ (j,5) _ (79) 

It is worth noticing that if the coefficient B were negative, then one should consider the decay rate associated to 
the process (j)(j> 4'4'4>j meaning that one would get a different temperature dependence for Fph, and thus, also for ^3. 

With the low energy effective field theory, one can as well compute the first corrections to the bulk viscosity 
coefficients for finite s-wave scattering length a. We have to evaluate how all the various ingredients needed in the 
computation are corrected in the presence of this scale breaking effect. 

The decay rate is modified as 

rCB^Fph('l-^jy)+0(y^), (80) 



where y has been defined in Eq. (|65p . 
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The expressions of Ii and I2 arc also affected and wc find that 



/f^ = h(l-'-^^^^^y]+Oiy^), (81) 



2ci + 3c2 

I?^ - ^T^(.^C(3)-7C(5))^!giy + 0(,^). (82) 
Upon substituting these expressions in Eqs. ((26|) . (pT]) and ((28|) we find that to leading order in y ^ 1/a 

Cr ^ -264.7 c. (c, + ^c.)^,, (83) 

7^3^3/2 

C2 ^ 19.0 cl^j^y\ (84) 
Mo' 

C3 . 369« (i 'f (ll±fcl!r3 (1 - 55£l±ip,) . ,85, 

V/io/ ™ V 8ci + 12c2 / 

Notice that the first non vanishing correction to is of the order of 1/a^, meaning that this bulk viscosity coefficient 
has to be neglected in the hydrodynamic equations ([5]), ([HI) and (jlO[) . where we are retaining terms of order 1/a only. 

VI. CONCLUSION 

We have derived the expressions for the three independent bulk viscosity coefficients of a non-relativistic Fermi 
superfiuid close to unitarity. In doing this we have first derived general expressions for these three transport coefficients 
assuming that the leading contribution comes from the superfiuid phonons, with dispersion law containing both a 
linear and a cubic term in momentum. Our computations are valid at very low temperatures, because for temperatures 
close to the superfiuid transition the contribution of other quasiparticles might also be important. 

In agreement with the outcome of Ref. [l^ we find by explicit calculation that at unitarity the bulk viscosity 
coefficients and ^2 vanish, whilst 

C,^ 3695Ah'(^ y^'^£L±M!(fc^^)^ (86) 

and we have restored here both the Planck and Boltzmann constants. In evaluating ^3 we have considered that the 
leadingcoUisional process is the phonon decay (jj — > </)(/>, which is kinematically allowed according to the results of 
Refs. (35l . [36| . If this were not the case, this transport coefficient would be dominated by the process (fxp — > 
resulting in a very different temperature dependence. 

Let us point out that in the computation of the shear viscosity for the same system [l^ , it was considered that the 
decay 4) ^ (fxj) was not kinematically allowed, resulting in a dependence r] oc l/T^. Considering this decay process 
for the computation of the shear viscosity, would change its temperature dependence into 77 oc l/T, as it occurs for 
^He (37| . Wc thus realize that it is very important to determine with good precision the phonon dispersion law. as 
the explicit values and temperature dependence of all the transport coefficients are extremely sensitive to its form. 

We have also evaluated corrections to the bulk viscosities due to conformal breaking terms, finding 

3 \e{kBTf 



Cr - -264.7;.C2^c, + -C2j^^,, (87) 
^ 19.0 , e'\ksTf . .... 

C2 - -Tr^2 — y : («») 

3695.4.^ r ^-^^^^(k.Tf (l - 'T!^]T:^) . (89) 

J .-2,3/2 

We have restricted the computation to the leading correction in the parameter y = h "^^'2^^ 1 which measures the 
departure from the conformal limit caused by a large but finite value of the scattering length a. 

The presence of non-vanishing bulk viscosity coefficients might be experimentally detectable. As an example it 
infiuences the propagation of first and second sound in a superfiuid [isj . The damping of first sound, ai, depends on 



yUo/ V 8ci + 12C2 
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the shear viscosity and on Ca-, whereas the damping of second sound, a2, depends on all the dissipativc cocfRcicnts. 
Since shear viscosity is much larger than the others it will give the leading contribution to the sound absorption 
coefRcients. However, one can show that in the unitary limit 

and therefore this combination of the absorption coefficients is independent of 77. Here is the frequency of the 
sound oscillation and k is the thermal conductivity. In order to evaluate the thermal conductivity one has to consider 
processes that change the total momentum of phonons, meaning that one has to consider the scattering of phonons 
with different (quasi)particles, see e.g. [ssj . A more detailed description of the thermal properties of unitary supcrfluids 
will be presented elsewhere. 

It might be interesting to explore the possibility of measuring the bulk viscosities in trapped Fermi superfluids 
through the study of dipole and breathing modes. The two-fluid equations of Landau predict the presence of different 
hydrodynamic modes, according to whether the superfluid and normal components oscillate in-phase, or move out- 
of-phase [S^. It has been suggested that the out-of-phase oscillations might be experimentally detected [i^l. If this 
is the case, its study might lead to the determination of the transport coefficients studied here, in the same way that 
the shear viscosity can be extracted studying the in-phase breathing mode [4l[. 
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